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We present the first measurement of the angle (j>3 of the unitarity triangle using a binned model- 
independent Dalitz plot analysis technique of B —> DK ± , D — > Kstt + tt~ decay chain. The method 
is based on the measurement of parameters related to the strong phase of D — s> Kgir + iT~ amplitude 
performed by the CLEO collaboration. The analysis uses full data set of 772 x 10 6 BB pairs 
collected by the Belle experiment at T(4S I ) resonance. We obtain 3 = (77.3±i|;| ± 4.2 ± 4.3)° and 
the suppressed amplitude ratio rs — 0.145 ± 0.030 ± 0.011 ± 0.011. Here the first error is statistical, 
the second is experimental systematic uncertainty, and the third is the error due to precision of 
strong phase parameters obtained by CLEO. This result is preliminary. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 

The angle 03 (also denoted as 7) is one of the least 
well-constrained parameters of the Unitarity Triangle. 
The theoretical uncertainties in 03 determination are ex- 
pected to be negligible, and the main difficulty in its 
measurement is a very low probability of the decays in- 
volved. 

The measurement that currently dominates the 03 sen- 
sitivity uses B^ — > DK^ decay with the neutral D me- 
son decaying to 3-body final state such as KgTr + n~ 
The weak phase 03 appears in the interference between 
b — > cus and b — > ucs transitions and can be measured in 
the Dalitz plot analysis of the D decay. This method re- 
quires the knowledge of the amplitude of D° —> A'g7r + 7r~ 
decay including its complex phase. The amplitude is ob- 
tained from the model that involves isobar and K-matrix 
descriptions of the decay dynamics, and thus results in 
the model uncertainty of the 03 measurement. In the 
latest model-dependent Dalitz plot analyses performed 
by BaBar and Belle, this uncertainty ranges from 3° to 

9° nm. 

The method to eliminate the model uncertainty us- 
ing the binned Dalitz plot analysis has been proposed by 
Giri et al. The information about the strong phase 
in D° — > Kgir + ir~ decay is obtained from the decays 
of quantum-correlated D° pairs produced in ^(3770) — > 
D°D° process. As a result, the model uncertainty is sub- 
stituted by the statistical error related to the precision 
of strong phase parameters. The method has been fur- 
ther developed in [t| [13] where the experimental feasi- 
bility of the method has been shown and the analysis 
procedure has been proposed to optimally use the avail- 
able B decays and correlated D° pairs. Recently, the 
measurement of strong phase in D° — > K s ir + Tr~ and 
D° — > K Q S K + K~ decays has been performed by CLEO 
collaboration [n], (T^. In this paper, we report the first 
measurement of 03 with the model-independent Dalitz 



plot analysis of D —> Kgir + ir~ decay from the mode 
_B ± — » DK^, based on a 711 fb _1 data sample (corre- 
sponding to 772 x 10 6 BB pairs) collected by the Belle 
detector at the KEKB asymmetric e + e~ factory using 
the CLEO measurement. These results are preliminary. 

II. THE MODEL-INDEPENDENT DALITZ 
PLOT ANALYSIS TECHNIQUE 

The amplitude of the B+ DK+, D K Q s ir+ir- 
decay is given by the interference of B + — > D°K + and 
B+ D°K+ amplitudes: 

A B =A + r B e l{&B+ ' l ' 3) A, (1) 

where A = A(m 2 Ksw+ ,m 2 Ksn _) = A(m 2 + ,m 2 _) is 

the amplitude of the D° — > Ag7r + 7r~ decay, A = 
A(m^_, m 2 _) is the amplitude of the D° — > A'g7r + 7r~ decay 
(A(m\,m 2 _) = A(m 2 _,m 2 t _) in the case of CP conserva- 
tion in D decay), r B is the ratio of the absolute values 
of the two interfering amplitudes, and S B is the strong 
phase difference between them. The Dalitz plot density 
of the D decay from B + — > DK + is given by 

Pb = \A B \ 2 =\A + r B e l{SB+ ^A\ 2 = 

r~= ( 2 ) 

P + r 2 B P + 2VPP(x + C + y+S), 

where 

x + = r B cos{5 B + 03 ) ; y + = r B sm(5 B + 03) ■ (3) 

The functions C = C{m 2 lr ,rn 2 _) and S = S(m\,m 2 _) 
are the cosine and sine of the strong phase difference 
8d = argA — &rgA between the D° —> Kgir + ir~ and 
D" A^tt+tt- amplitude^]: 

C = cos 5 D{rn 2 + ,m 2 _)-, S = sm8 D (m 2 + ,m 2 _) . (4) 



* deceased 



This paper follows the convention for strong phases in D decay 
amplitudes introduced in Ref. ficl . 
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The equations for the charge-conjugate mode B~ — > 
DK~ arc obtained with the substitution 3 — > — 3 and 
Ah4; the corresponding parameters of the admixture 
of the suppressed amplitude are 

x- = r B cos(S B - fa) 5 V- = Tb sm(5 B - <fo) ■ (5) 

Using both B charges, one can obtain 03 and 5 B sepa- 
rately. 

Up to this point, the description of the model- 
dependent and model-independent techniques is the 
same. The model-dependent analysis deals directly with 
the Dalitz plot density, and the functions C and S 
are obtained from model assumptions in the fit of the 
D° — > KgTr + ir~ amplitude. In the model-independent 
approach, the Dalitz plot is divided into 27V bins sym- 
metric under the exchange m 2 _ o m 2 + . The bin index 
"i" ranges from —TV to TV (excluding 0); the exchange 
m\ <-> m 2 _ corresponds to the exchange i -H* — i. The 
expected number of events in the bin "z" of the Dalitz 
plot of D from B+ -> DK+ is 



Ki 



rlK. 



2^K l K^ i (x + c i + y + Si) 



(6) 

where h B is a normalization constant and Ki is the num- 
ber of events in the z-th bin of the Dalitz plot of the D 
meson in a flavor eigenstate (obtained using D** — > Dtt 
sample). The terms c; and Si include information about 
the cosine and sine of the phase difference averaged over 
the bin region: 



/ \A\\A\ cos S D dV 
I \A\ 2 dV j \A\ 2 dV 



(7) 



Here V represents the Dalitz plot phase space and T>i is 
the bin region over which the integration is performed. 
The terms are defined similarly with cosine substituted 
by sine. 

The absence of CP violation in D decay requires a = 
c-i and Si = — s_». The values of a and s, terms can 
be measured in the quantum correlations of D pairs by 
charm-factory experiments operated at the threshold of 
DD pair production [ill ■ The wave function of the 
two mesons is antisymmetric, thus the four-dimensional 
density of two correlated D — > Kgir + ir~ Dalitz plots is 

\A COII (m 2 + ,m 2 ^,m' 2 ,m' 2 )\ 2 = \A{A 2 -^A 2 \ 2 = 



P1P2 + P1P2 - 2\l P l P 2 P 1 P 2 {C x C 2 + SxS 2 ) , 

(8) 

where the indices "1" and "2" correspond to the two 
decaying D mesons. In the case of a binned analysis, the 
number of events in the region of the (K^ir + 'K~) 2 phase 
space described by the indices "i" and "j" is 



Ma =KiK-j + K-iKj 



Another possibility is to use decays of D in a CP- 
cigenstate to KgTr + n~ tagged by the other D decaying to 
a CP-odd or CP-even final state. This allows to obtain 
constraints on the value of c;. 

Once the values of the terms c, and s, are known, the 
system of equations (j6|) contains only three free parame- 
ters (x, y, and h B ) for each B charge, and can be solved 
using a maximum likelihood method to extract the value 
of 3 . 

We have neglected charm mixing effects in decays of 
D from both the B^ — > DK^ process and in quantum- 
correlated DD production. It has been shown [l3| that 
although charm mixing correction is of the first order in 
the mixing parameters x,y, it is numerically small (of 
the order 0.2° for x,y ~ 0.01) and can be neglected at 
the current level of precision. The future precision mea- 
surements of 4>3 can account for both charm mixing and 
CP violation (both in mixing and decay) once the corre- 
sponding parameters are measured. 

Note that technically the system ([6]) can be solved 
without external constraints on c; and s, for TV > 2. 
However, due to the small value of r B , there is very little 
sensitivity to the Cj and Sj parameters in B ± — s- DK^ 
decays, which results in a reduction in the precision on 
03 that can be obtained S. 



III. CLEO INPUT 

The procedure of binned Dalitz plot analysis should 
give the correct results for any binning. However the 
statistical accuracy depends strongly on the amplitude 
behavior across the bins. Significant variations of the 
amplitude within a bin results in loss of coherence in the 
interference term. This effect becomes especially signifi- 
cant with limited statistics when a small number of bins 
has to be used. Better statistical precision is reached for 
the binning with the phase difference between the D° 
and D° amplitudes varying as little as possible For 
the optimal precision, one also has to take the variations 
of the absolute value of the amplitude into account and 
the contribution of background events. The procedure to 
optimize the binning from the point of view of 03 statisti- 
cal precision has been proposed i n fioll and generalized to 
the case with the background in [lj]. It has been shown 
that as little as 16 bins are enough to reach the statistical 
precision just 10-20% smaller than in the unbinned case. 

The optimization of binning sensitivity uses the am- 
plitude of D — > Kgir + ir~ decay. It should be noted, 
however, that although the choice of binning is model- 
dependent, a poor choice of model results only in the loss 
of precision but not in the bias of measured parameters. 
CLEO measured Ci , Si terms for four different binnings 



with 16 bins (bin index i £ (—8, . 



.8)): 



-2 V /K l K_iKj K_j Cj + SiSj). 



(9) 



1. The binning optimized for statistical precision ac- 
cording to the procedure from [l(| (see Fig.[T]). The 
effect of the background is not taken into account. 



5 




FIG. 1. Optimal binning of D K%-k + k~ Dalitz plot. 



The amplitude is taken from the BaBar measure- 
ment 

2. Same as above, but optimized for the analysis with 
high background in B data (e. g. at LHCb) 

3. The binning with bins equally distributed in the 
phase difference A5d between the D° and D de- 
cay amplitudes, with the amplitude from the BaBar 
measurement [4|. 

4. Same as above, but with the amplitude from the 
Belle analysis 

Our analysis uses the optimal binning shown in Fig. Q] 
(option 1) as the baseline since it offers better statistical 
accuracy. In addition, we use the equal phase difference 
binning (A^p-binning, option 3) as a cross-check. 

The results of the CLEO measurement of Cj,Sj terms 
for the optimal binning are presented in Table |TJ The 
same results in graphical form arc shown in Fig. [2] 
The comparison with the a, Si calculated from the Belle 
model is presented, and shows a reasonable agreement 
between the Belle model and measurement. 



IV. ANALYSIS PROCEDURE 

The key equation of the analysis (JSJ holds in the ideal 
situation with no background, constant efficiency (inde- 
pendent of Dalitz plot variables), and no cross- feed be- 
tween the bins due to momentum resolution and radiative 
corrections. In this section, we outline the procedure of 
the analysis taking the above mentioned effects into ac- 
count. 



TABLE I. Values of the phase terms d,Si for the optimal 
binning measured by CLEO [12], and calculated from Belle 
D — > Kg-K + ir~ amplitude model. The 'x 2 /ndf of the agree- 
ment between the measured and predicted ct,Si is 18.6/16. 





CLEO measurement 


Belle model 


Cl 


-0.009 ±0.088 ±0.094 


-0.039 


C2 


±0.900 ±0.106 ±0.082 


±0.771 


C:>, 


±0.292 ±0.168 ±0.139 


±0.242 


C4 


-0.890 ±0.041 ±0.044 


-0.867 


Cg 


-0.208 ±0.085 ±0.080 


-0.246 


C6 


±0.258 ±0.155 ±0.108 


±0.023 


C7 


±0.869 ±0.034 ±0.033 


±0.851 


C8 


±0.798 ±0.070 ±0.047 


±0.662 


Si 


-0.438 ±0.184 ±0.045 


-0.706 


Si 


-0.490 ±0.295 ±0.261 


±0.124 


S3 


-1.243 ±0.341 ±0.123 


-0.687 


Si 


-0.119 ±0.141 ±0.038 


-0.108 


S 5 


±0.853 ±0.123 ±0.035 


±0.851 


S6 


±0.984 ±0.357 ±0.165 


±0.930 


S7 


-0.041 ±0.132 ±0.034 


±0.169 


s 8 


-0.107 ±0.240 ±0.080 


-0.596 
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FIG. 2. Comparison of phase terms Ci,Si for the optimal 
binning measured by CLEO, and calculated from the Belle 
D — > K%n + n~ amplitude model. 



A. Efficiency profile 

The effect of non- uniform efficiency over the Dalitz plot 
is canceled out by using the flavor-tagged D sample with 
the same kinematic properties as the sample for the signal 
B decay. This approach allows of removal of the system- 
atic error associated with the possible inaccuracy in the 
description of the detector acceptance in the Monte- Carlo 
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(MC) simulation. 

We note that the Equations © and ((5J do not change 
after the transformation p — > ep, if the efficiency profile 
e(m^_,m?_) is symmetric: e(m^_,m 2 _) = e(m^,m^_). It 
means that if the efficiency profile is the same in all three 
measurements involved (flavor D, correlated "0(3770) — > 
DD, and D from B — > DK), the resulting measurement 
will be unbiased even if no efficiency correction is applied. 

We match the Dalitz plot efficiency profiles for the fla- 
vor D to the one from B — > DK by taking the flavor- 
tagged D mesons with the same average momentum as 
the D mesons from B — > DK. The center-of-mass (CM) 
D momentum distribution for B — > DK decays is prac- 
tically uniform in the narrow range 2.1 GeV/c < pd < 
2.45 GeV/c. We assume that the efficiency profile de- 
pends mostly on the D momentum and take the flavor- 
tagged sample with the average momentum of po = 2.3 
GeV/c (we use a wider range of D momenta than in 
B — > DK to increase the statistics). The assumption 
that the efficiency profile depends only on the D momen- 
tum is tested using MC simulation, and the remaining 
difference is treated as the systematic uncertainty. 

While calculating Ci,Si, CLEO applies an efficiency 
correction, therefore the values reported in their anal- 
ysis correspond to flat efficiency profile. To use the Cj, Si 
values in the fa analysis, they have to be corrected for 
the Belle efficiency profile. This correction cannot be 
performed in a completely model-independent way, since 
the correction terms include the phase variation inside 
the bin. Fortunately, the calculations using the Belle 
D —> KgTr + Tr^ model show that this correction is negli- 
gible even for very large non-uniformity of the efficiency 
profile. The difference between the uncorrected ci,Si 
terms and those corrected for the efficiency, calculated 
using the efficiency profile parametrization used in the 
605 fb _1 analysis Q, does not exceed 0.01, i. e. it is 
negligible compared to the statistical error. 

B. Momentum resolution 

Finite momentum resolution leads to migration of 
events between the bins. In the binned approach, this 
effect can be corrected for in a non-parametric way. The 
migration can be described by the linear transformation 
of the number of events in bins: 

Ni = J2 a ikN k , (10) 

where Nk is the number of events the bin k would con- 
tain without the cross-feed, and N[ is the reconstructed 
number of events in the bin i. The cross-feed matrix oak 
is nearly unit: 1 for i ^ k. It is obtained from 

the signal MC simulation with the amplitude measured in 
Belle's 605 ftr 1 analysis Q. In the case of D ->■ K%ir + -K~ 
decay from B, the cross- feed depends on the parameters 
x,y. We assume that this effect is small and neglect it. 

Migration of events between the bins occurs also due 
to final state radiation (FSR). The Ci,Si terms in the 



CLEO measurement are not corrected to FSR; we there- 
fore do not simulate FSR to obtain the cross-feed matrix 
to minimize the bias due to this effect. Comparison of the 
cross-feed with and without FSR shows that this effect 
is negligible. 



C. Fit procedure 

The background contribution has to be accounted for 
in the calculation of the values Ni and K{. Statistically 
the most effective way of calculating the number of signal 
events (especially in the case of Ni, where the statistics 
is a limiting factor) is to perform the unbinned fit in the 
event selection variables for the events in each bin "i" of 
the Dalitz plot. 

Two different approaches are used in this analysis. In 
the first one, we fit the data distribution in each bin sep- 
arately, with the number of events for signal and back- 
grounds as free parameters. Once the numbers of events 
in bins Ni are found, we use them in Eq. |5]to obtain the 
parameters x±,y±. Technically it is done by minimizing 
the negative logarithmic likelihood of the form 

-21og£(.T,y) = -2j2^gp({Ni)(x,y),Ni,a Ni ), (11) 

i 

where (N)(x,y) is the expected number of events in the 
bin i obtained from Eq. [6l N and cr^ are the observed 
number of events and its error obtained from the data fit. 
If the probability density function (PDF) p is Gaussian, 
this procedure translates to the x 2 fit. 

The procedure described above does not make any as- 
sumptions on the Dalitz distribution of the background 
events, since the fits in each bin are independent. Thus 
there is no associated systematic uncertainty. However, 
in the case of low number of events and many background 
components this can be a limiting factor. Another solu- 
tion is to use the combined fit with a common likelihood 
for all bins. Relative numbers of background events in 
bins in such a fit can be constrained externally from e. g. 
MC sample. In addition, in the case of the combined fit , 
the two-step procedure of first extracting the numbers 
of signal events, and then using them to obtain (x, y) is 
not needed — the expected numbers of events (Ni) as 
functions of (x, y) can be plugged directly into the likeli- 
hood. Thus the variables (x, y) become free parameters 
of the combined likelihood fit, and the assumption of the 
Gaussian distribution of the number of signal events is 
not needed. 

Both approaches (separate fits in bins, and the com- 
bined fit) are tested with the control sample and the MC 
simulation. We choose the combined fit approach as the 
baseline, but the approach with separate fits in bins is 
also used: it allows to clearly demonstrate the CP asym- 
metry of the number of events in bins. 
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V. EVENT SELECTION 

We use a data sample of 772 x 10 6 BB pairs collected 
by the Belle detector. The decay chains B + — > DK + and 
are selected for the analysis. The neutral I? 
meson is reconstructed in the Kgir + ir~ final state in all 
cases. We also select decays of D* + — > D Q n + produced 
via the e + e~ — > cc continuum process as a high-statistics 
sample to determine the terms related to flavor-tagged 
L>° -> Xg7r+7T- decay. 

The Belle detector is described in detail elsewhere 
fbH [l5j . It is a large- solid- angle magnetic spectrometer 
consisting of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC) for charged particle tracking 
and specific ionization measurement (dE/dx), an array 
of aerogel threshold Cherenkov counters (ACC), time- 
of- flight scintillation counters (TOF), and an array of 
CsI(Tl) crystals for electromagnetic calorimctry (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux return located 
outside the coil is instrumented to detect Kl mesons and 
identify muons (KLM). 

Charged tracks are required to satisfy criteria based 
on the quality of the track fit and the distance from the 
interaction point (IP). We require each track to have a 
transverse momentum greater than 100 MeV/c, and the 
impact parameter relative to the IP of the beams less 
than 2 mm in the transverse and less than 10 mm in 
longitudinal projections. Separation of kaons and pi- 
ons is accomplished by combining the responses of the 
ACC and the TOF with the dE/dx measurement from 
the CDC. Neutral kaons are reconstructed from pairs of 
oppositely charged tracks with an invariant mass M W7r 
within 7 MeV/c 2 of the nominal Kg mass, flight dis- 
tance from the IP in the plane transverse to the beam 
axis grater than 0.1 mm, and the cosine of the angle 
between the projections of flight direction and its mo- 
mentum greater than 0.95. 

To determine the terms for flavor-tagged D° — > 
K s tt + it~ decay we use mesons produced via the 

e + e _ — > cc continuum process. The flavor of the neu- 
tral D meson is tagged by the charge of the slow pion 
in the decay D* + —> D°tt + . The slow pion track is re- 
quired to originate from the D° decay vertex to improve 
the momentum and angular resolution. The selection of 
signal candidates is based on two variables, the invariant 
mass of the neutral D candidates Mb = M K o 7T+7T ^- and 
the difference of the invariant masses of the D* ± and the 
neutral D candidates AM = Af ( - i ^o 7r + 7r -- )£)7r — M K o n + n -. 
We retain the events satisfying the following criteria: 
1800 MeV/c 2 < M fl < 1920 MeV/c 2 and AM < 0.15 
MeV/c 2 . We also require the momentum of the D° can- 
didate in the CM frame pn to be greater than 1.5 GeV/c. 
About 15% of selected events contain more than one D*^ 
candidate that satisfies the requirements above; in that 
case we keep only one randomly selected candidate. 

Selection of B^ —> DK^ and i? ± — > Dir^ samples 
is based on the CM energy difference AE = Ei — 



£ bcam an(1 the bea m-constrained B meson mass Mb c = 
V E b ca .m -(Y.Pi) 2 -, where E hcam is the CM beam energy, 
and Ei and pi are the CM energies and momenta of the B 
candidate decay products. We select events with Mb c > 
5.2 GeV/c 2 and \AE\ < 0.18 GeV for further analysis. 
We also impose a requirement on the invariant mass of 
the neutral D candidate \M K o^+^- -M D o | < 11 MeV/c 2 . 

Further separation of the background from e + e~ — > qq 
(q = u, d, s, c) continuum events is done by calculating 
two variables that characterize the event shape. One is 
the cosine of the thrust angle cosflthr, where #thr is the 
angle between the thrust axis of the B candidate daugh- 
ters and that of the rest of the event, calculated in the 
CM frame. The other is a Fisher discriminant F com- 
posed of 11 parameters the production angle of the 
B candidate, the angle of the B thrust axis relative to 
the beam axis, and nine parameters representing the mo- 
mentum flow in the event relative to the B thrust axis in 
the CM frame. 

In both flavor D° and B ± -> DK ± (B ± -> Dtt ± ) sam- 
ples the momenta of the tracks forming a D° candidate 
are fitted to the nominal D° mass in the calculation of 
the Dalitz plot variables. 

VI. FLAVOR- TAGGED SAMPLE D*+ ->■ £>%+, 
D -> K° s -r + -r- 

The numbers of events Ki in bins of the flavor- 
tagged D — > Kgir + ir~ decay are obtained from the two- 
dimensional unbinned fit of the distribution of variables 
Ma and AM. The fits in each Dalitz plot bin are done 
independently. We take the candidates in the CM D 
momentum range 1.8 GeV/c < pn < 2.8 GeV/c. It 
provides the same average po as in B — > DK decays 
(po = 2.3 GeV/c) to reduce the influence of the effi- 
ciency profile on 03 measurement (see Section TlV Ap . 

The fit uses the signal PDF and two background com- 
ponents: purely random combinatorial background and 
the background with real D° and random slow pion track. 
The signal distribution is a product of the PDF's for Mjj 
(triple Gaussian) and AM (sum of bifurcated Student 
distribution and bifurcated Gaussian distribution). The 
combinatorial background is parametrized by the linear 
function in Mp and by the function with a kinematic 
threshold at the ir + mass in AM. The random slow 
pion background is parametrized as a product of signal 
Md distribution and the combinatorial AM background 
shape. 

The parameters of the signal and background distribu- 
tions are obtained from the data fit. The parameters of 
the signal PDF are constrained to be the same in all bins. 
The free parameters in each bin are the number of signal 
events Ki , the parameters of the background distribution 
and fractions of the background components. 

The fit results of the flavor sample for the whole Dalitz 
plot are shown in Fig. [3] The number of signal events 
calculated from the integral of the signal distribution is 
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FIG. 3. Projections of the flavor-tagged D* + Z>°7r + , D -¥ 
Ksiv + tt~ data with 1.8 GeV/c < p D < 2.8 GeV/c onto (a) 
Mb and (b) AAf variables. Histograms show fitted signal and 
background contributions, points with the error bars are the 
data. Full D — ¥ Kgii + -K~ Dalitz plot is used. 



TABLE II. Numbers of events in Dalitz plot bins for the 
flavor-tagged D* + — > D°ir + , D — > Kg-K + iv~ sample with 
1.8 GeV/c <p D < 2.8 GeV/c. Results of the 2D (AM, M D ) 
fit to data. 



Bin i 


Kiipn G (1.8,2.8)) 


-8 


26450 ± 181 


-7 


22476 ± 196 


-6 


1765 ± 68 


-5 


13146 ± 143 


-4 


26482 ± 202 


-3 


1601 ± 58 


-2 


1827 ± 63 


-1 


8770 ± 124 


1 


43261 ± 255 


2 


58005 ± 268 


3 


62808 ± 274 


4 


44513 ± 253 


5 


21886 ± 177 


6 


28876 ± 197 


7 


48001 ± 265 


8 


9279 ± 125 


Total 


426938 ± 825 



426938 ± 825, the background fraction in the signal box 
\M D -m D o\ < 11 MeV/c 2 , 144.5 MeV/c 2 < AM < 146.5 
MeV/c 2 is 10.1 ± 0.1%. The numbers of events in bins 
are shown in Table |TTJ 



VII. SELECTION OF B ± -> Dvr* AND B ± -> DK ± 
SAMPLES 

The decays B ± — > DK^ and B ± — > Dtt^ have similar 
topology and background sources and their selection is 
performed in a similar way. The mode B ± — > Dir^ has 



an order of magnitude larger branching ratio and small 
amplitude ratio rg ~ 0.01 due to small ratio of the weak 
coefficients \V ub V* d \ / \V cb V* d \ ~ 0.02 and additional color 
suppression factor as in the case of B ± — > DK^. This 
mode is used as a control sample to test the procedures 
of the background extraction and Dalitz plot fit. Also, 
the resolution scale factors and Dalitz plot structure of 
some background components are constrained from the 
control sample and used in the signal fit. 

Extraction of the number of signal events is performed 
by fitting the 4D distribution of variables Mb c , AE, 
cos 6> t hr and J ' . The fit to B ± — > sample uses three 
background components in addition to the signal PDF. 
These are: 



• Combinatorial background from 
cess, where q = (u, d, s, c). 



qq pro- 



• Random BB background, where the tracks form- 
ing the B ± — > Dtt^ 1 candidate come from decays 
of both B mesons in the event. The number of 
possible B decay combinations that contribute to 
this background is large, therefore both the Dalitz 
distribution and (Mb c , AE) distribution are quite 
smooth. 

• Peaking BB background, where all tracks forming 
the B^ — > D7r ± candidate come from the same B 
meson. This kind of background is dominated by 
the B —> D*tt decays with lost ir or 7 from D* 
decay 

The i? ± — > DK^ fit includes in addition the background 
from B^- — > Dn^ decays with pion misidcntified as kaon. 

The PDF for the signal parametrization (as well as 
for each of the background components) is a product of 
(M hc ,AE) and (cos t hr, J 7 ) PDFs. The (M hc ,AE) PDF 
is a sum of two 2D Gaussian distributions (core and wide) 
with the correlation between Mb c and AE. We use com- 
mon parametrization for the cos^thr,^ 7 distribution for 
signal and all background components. The distribution 
is parametrized by the sum of two functions (with differ- 
ent coefficients) of the form 



p(x,F) =exp(Cia; + C 2 x 2 + C 3 x 3 )x 

G(T, F (x), cr FL (x), a FR (x)), 



(12) 



where x = cos(9 t hr, G(x,F,or,oii) is the bifurcated 
Gaussian distribution with the mean F and the widths 
<7l and an, and functions Fq, <tfl and <jfr are the poly- 
nomials which contain only even powers of x. 



Combinatorial background from continuum e + e~ — > 
qq production is obtained from the experimental sam- 
ple collected at the CM energy below T(45) resonance 
(off- resonance data). The parametrization in variables 
(cos #thr, -T 7 ) is the same as described for the signal PDF. 
The parametrization in (Mbc AE) is the product of ex- 
ponential distribution in AE and the empirical shape 
proposed by the Argus collaboration [l7j in Mb c : 

Pcomb(M hc , AE) = exp(-aAE)M hcy /yexp(-cy), (13) 
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where y = 1 — Mb c /i?boam, and -Ebeam is the CM beam 
energy. 

Random BB background is obtained from generic MC 
sample. Generator information is used to select only the 
events where the candidate is formed from tracks com- 
ing from both B mesons. The (Mb c , AE) distribution 
of this background is parametrized by the sum of three 
components: 

• product of exponential (in AE) and Argus (in Mb c ) 
functions, as for continuum background, 

• product of exponential in AE and bifurcated Gaus- 
sian distribution in M\, c , where the mean of the 
Gaussian distribution is linear as a function of AE. 

• two-dimensional Gaussian distribution in AE, Mbc 
with correlation (asymmetric in M^ c ). 

Peaking BB backgrounds are parametrized by the 
same function as the random BB background. The 
background coming from B + B~ and B°B° decays is 
treated separately in (Mbc, AE) variables, while the com- 
mon (cos (9thr, E) distribution is used. In the case of 
B ± -> DK± fit, we take the B ± -> Dn ± events with 
pion misidentificd as kaon as a separate background cat- 
egory. The distributions of Mb c , AE and cos# t hr, F vari- 
ables are parametrized in the same way as for the signal 
events and are obtained from MC simulation. 

The Dalitz plot distributions of the background com- 
ponents are discussed in the next section. Note that the 
Dalitz distribution is described by the relative number 
of events in bins. The numbers of events in bins can 
be taken as free parameters in the fit, thus there will 
be no uncertainty due to Dalitz plot description of the 
background in such an approach. This procedure is jus- 
tified for the background that is well separated from the 
signal (such as peaking BB background in the case of 
B ± — > Z?7r ± ), or if the background is constrained by 
much larger number of events than the signal (such as 
the continuum background). 

The results of the fit to B ± -> Dn ± and B ± -> DK ± 
data with the full Dalitz plot taken are shown in Figs 0] 
and [5j respectively. We obtain a total of 19106 ± 147 
signal B ± — > Dir* events and 1176 ± 43 signal B ± — > 
DK events — 55% more than in the 605 fb _1 model- 
dependent analysis Q . The improvement partially comes 
from larger integrated luminosity, and partially from the 
better selection efficiency due to improved tracking pro- 
cedures. 



VIII. DATA FITS IN BINS 

The data fits in bins for both B ± -> Z?tt ± and B ± -> 
DK ± are performed with two different procedures: sep- 
arate fits for the number of events in bins, and the com- 
bined fit with the free parameters (x,y), as discussed in 
Section HV CI The combined fit is used to obtain the final 
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FIG. 4. Projections of the B ± — > Di* data onto (a) M bc , 
(b) AE, (c) cos #thr and (d) T variables. Histograms show 
fitted signal and background contributions, points with the 
error bars are the data. Full D — ► KgTv + ir~ Dalitz plot is 
used. 
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values for (x, y), while the separate fits provide the cross- 
check of the fit procedure and a way to visualize the CP- 
violating effect. The study with MC pseudo-experiments 
is performed to check that the observed difference in the 
fit results between the two approaches agrees with the 
expectation. 

In the case of separate fits in bins, we first perform 
the fit to all events in the Dalitz plot. The fit uses back- 
ground shapes fixed from the generic MC simulation of 
continuum and BB decays. The signal shape is fixed 
from the signal MC sample except that we float the mean 
values of AE and M^ c and width scale factors. As a next 
step, we fit the 4D (Mb c , AE, cos# t hr: F) distributions in 
each bin separately. The free parameters of each fit are 
the number of signal events, and the number of events in 
each background category. 

The numbers of signal events in bins extracted from 
the fits are given in Table Mil These numbers are used in 
the fit to extract (x,y) using © after the cross-feed and 
efficiency correction for both N{ and JQ. Figure [6] illus- 
trates the results of this fit. The numbers of signal events 
in bins separately for B + and B~ are shown in Fig. Eta) 
together with the numbers of events in the flavor sample 
(appropriately scaled). The difference in the number of 
signal events shown in Fig. IH^b) does not reveal CP vi- 
olation. Figures [S^c) and (d) show the difference of the 
numbers of signal events for B + (B~) data and scaled 
flavor sample, both for the data and after the (x, y) fit. 
The x 2 /ndf is reasonable for both the (x, y) fit and for 
the agreement with the purely flavor-specific amplitude. 

Unlike B ± -> Dn ± , the B ± -> DI^ sample shows 
significantly different numbers of events in bins of B + 
and B~ data (see Fig. \T[b) and Table [[V). The prob- 
ability to obtain this difference as a result of statistical 
fluctuation is 0.42%. This number can be taken as the 
model-independent measure of the CP violation signifi- 
cance. The significance of fa being nonzero is in general 
smaller since fa ^ results in a specific pattern of charge 
asymmetry. The fit of the numbers of events to the ex- 
pected pattern described by the parameters (x, y) shows 
a good quality[7tc,d), i. e. is consistent with the hypoth- 
esis that the observed CP violation is solely explained by 
the mechanism involving nonzero fa. 

The default combined fit uses the constraint of the 
random BB background in bins from the generic MC, 
and x±,y± variables as free parameters. Fits to B + 
and B~ data are performed separately. Additional free 
parameters are the numbers of continuum and peaking 
BB backgrounds in each bin, fraction of the random 
BB background, and means and scale factors of the sig- 
nal Mbc and AE distributions. The values of x, y are 
then corrected for the fit bias obtained from MC pseudo- 
experiments. The value of the bias depends on the initial 
x and y values and is of the order 5x 10 -3 for B^ — > DK^ 
sample and less than 10 -3 for B^ — > Dtt^ sample. 



TABLE III. Numbers of events in Dalitz plot bins for the 
B — > Dtt , D —¥ K%-k + tt~ sample with the optimal bin- 
ning. Results of the independent 4D fits with variables 
(Mbc, AE, cos # t hr, F) fit to data. 



Bin i 






-8 


564.2 ± 25.3 


587.0 ± 25.7 


-7 


462.3 ± 23.8 


462.8 ± 23.9 


-6 


47.9 ± 7.7 


39.2 ± 7.2 


-5 


314.1 ± 19.0 


286.2 ± 18.2 


-4 


592.6 ± 26.5 


645.7 ± 27.8 


-3 


22.2 ±6.2 


27.2 ± 6.3 


-2 


42.7 ±7.6 


54.0 ± 8.7 


-1 


190.8 ± 15.4 


210.8 ± 16.3 


1 


959.2 ± 32.6 


980.2 ± 33.1 


2 


1288.7 ± 37.0 


1295.9 ± 37.1 


3 


1395.8 ± 38.4 


1352.2 ± 37.9 


4 


1045.5 ± 34.7 


1065.1 ±34.9 


5 


479.3 ± 23.3 


532.2 ± 24.5 


6 


623.7 ±26.0 


663.5 ± 26.7 


7 


1081.0 ±35.3 


1049.2 ± 34.8 


8 


210.0 ± 16.1 


212.1 ± 16.3 


Total 


9467.1 ± 103.6 


9639.1 ± 104.7 
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FIG. 6. Results of the fit of B — > Di ± control sample, 
(a) Numbers of events in bins of D — > A'g7r + 7r~ Dalitz plot: 
from B~ — > Dn~ (red), B + — > D% + (blue) and flavor sam- 
ple (histogram), (b) Difference of the number of events from 
B + — > Dn + and B~ —5- Dn~~ decays, (c) Difference of the 
number of events from B~ — > Dn~ and flavor sample (nor- 
malized to the total number of B~ —¥ Dtt~ decays): data 
(points with the error bars), and as a result of the (x,y) fit 
(horizontal bars), (d) Same for B + —¥ Dir + data. 



The values of x, y parameters and their statistical cor- 
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TABLE IV. Numbers of events in Dalitz plot bins for the 
B ± DK ± , D ->■ Kg 7T + 7r sample with the optimal 
binning. Results of the independent 4D fits with variables 
(M bc , AE, cos 6 tbv ,T) fit to data. 



Bin i 




N? 


-8 


49.8 ± 8.2 


37.8 ± 7.5 


-7 


42.2 ± 8.6 


24.9 ± 7.2 


-6 


0.0 ±1.9 


3.4 ±2.9 


-5 


9.6 ±4.5 


23.6 ± 6.2 


-4 


32.9 ±7.5 


42.1 ± 8.3 


-3 


3.5 ±2.8 


0.7 ±2.5 


-2 


11.3 ±4.1 


0.0 ± 1.3 


-1 


16.6 ± 5.4 


7.7 ±4.4 


1 


37.6 ±8.0 


65.1 ±9.9 


2 


68.6 ±9.6 


75.5 ±9.8 


3 


83.4 ± 10.1 


82.4 ± 10.2 


4 


49.3 ± 9.1 


86.5 ± 11.4 


5 


34.0 ± 7.3 


38.3 ± 7.6 


6 


34.8 ± 6.8 


41.9 ±7.5 


7 


70.8 ± 10.6 


46.4 ± 9.0 


8 


9.4 ±4.3 


14.2 ± 5.1 


Total 


574.9 ± 29.9 


601.6 ±30.8 
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FIG. 7. Results of the fit of B -> control sample, 

(a) Numbers of events in bins of D — >• _Klj7r + 7r~ Dalitz plot: 
from B~ — )• M" (red), B + — > DK + (blue) and flavor sam- 
ple (histogram), (b) Difference of the number of events from 
B + — > DK + and B~ — > DK~ decays, (c) Difference of the 
number of events from B~ — > DK~ and flavor sample (nor- 
malized to the total number of B~ — ¥ DK~ decays): data 
(points with the error bars), and as a result of the (x,y) fit 
(horizontal bars), (d) Same for B + — > DK + data. 
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FIG. 8. One-, two-, and three standard deviations levels for 



x,y fit of B -> DK mode. 



relations obtained from the combined fit are as follows: 

= -0.0045 ± 0.0087 ± 0.0050 ± 0.0026, 
y- = -0.0231 ± 0.0107 ± 0.0050 ± 0.0065, 

corr(x_,y_) = -0.189, 

x + = -0.0172 ± 0.0089 ± 0.0060 ± 0.0026, 
y+ = ±0.0129 ± 0.0103 ± 0.0060 ± 0.0065, 

corr(x+,y+) = -0.205 

(14) 

for i? ± — > £>7r ± control sample and 



X- = 


±0.095 ± 0.045 ± 0.014 ± 0.017, 


V- = 


+0.137±g;gp ± 0.019 ± 0.029, 


corr(x_,y_) = 


-0.315, 


x+ = 


-0.110 ± 0.043 ± 0.014 ± 0.016, 


y+ = 


-0.050i^5 ± 0.011 ± 0.021, 


corr(x+,y + ) = 


±0.059 



(15) 



for S ± — > DK ± sample. Here the first error is statisti- 
cal, the second error is the systematic uncertainty, and 
the third error is the uncertainty due to the errors of 
Ci,Si terms. The measured values of (x±,y±) with their 
likelihood contours arc shown in Fig. [8] 



IX. SYSTEMATIC ERRORS 

Systematic errors in the x, y fit are obtained for the 
default procedure of the combined fit with the optimal 
binning. The systematic errors are summarized in Ta- 
ble El 

The uncertainty of the signal shape used in the fit in- 
cludes the following sources: 
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• Choice of paramctrization used to model the shape. 
The corresponding uncertainty is estimated by us- 
ing the non-parametric (Keys) function instead of 
parametrized distribution. 

• Possible correlation between the (Mb c , A£) and 
(cos (9 t hr, E) distributions. To estimate its effect we 
use 4D binned histogram to describe the distribu- 
tion. 

• MC description of the (cos#thr, J~) distribution. Its 
effect is estimated by floating the parameters of the 
distribution in the fit to — > Dtt^ control sam- 
ple. 

• Dependence of the signal width on the Dalitz plot 
bin. The uncertainty due to this effect is esti- 
mated by performing the B^ — > Dir^ fit with the 
shape parameters floated separately for each bin, 
and then using the results in the fit to B^ — > DK ± 
data. 

We do not assign the uncertainty due to the difference 
in (Mbc, AE) shape between the MC and data since the 
width of the signal distribution is calibrated on B ± — > 
Dn^ data. 

In the uncertainty of the continuum background shape, 
the same four sources are considered as for the sig- 
nal distribution. The uncertainty due to the choice of 
parametrization is estimated similarly by using the Keys 
function. The effect of possible correlation between the 
(Mbc, AE) and (cos $thr, J 7 ) distributions is estimated by 
using the distribution split into the sum of two compo- 
nents (it, d, s and charm contributions) with independent 
(Mbc AE) and (cos t hr , J~) shapes. The uncertainty due 
to MC description of the (Mbc AE) and (cos#thr, J~) dis- 
tributions is estimated by floating their parameters in the 
B ± — > Dtt fit. To estimate the effect of possible corre- 
lation of the shape with the Dalitz plot variables we fit 
the shapes separately in each Dalitz plot bin. 

The uncertainties of the shapes of random and peak- 
ing BB backgrounds are estimated conservatively by per- 
forming the fit with AE > —0.1 GeV — this requirement 
rejects the peaking BB background and a large part of 
the random BB background. 

In the case of the fit of i? ± — > DK^ sample, 
the uncertainty of the B ± — > D7r ± background shape 
in (cos #thr, -E) variables is estimated by taking the 
(cos^thr,-? 7 ) shape for signal events. The Dalitz plot 
uncertainty is estimated by using the number of flavor- 
tagged events in bins (rather than the number of B — > 
Z?7r ± events used in the default fit). Uncertainties due 
to possible correlations are treated as in the case of the 
signal distribution. 

The uncertainty due to Dalitz plot efficiency shape ap- 
pears because of a difference in average efficiency for the 
flavor and B — > DK ± samples. The maximum differ- 
ence of 1.5% is obtained in the MC study. The uncer- 
tainty is obtained from the maximum of two quantities: 



• RMS of x and y from smearing the numbers of 
events in the flavor sample Ki by 1.5%. 

• Bias of x and y between the fits with and with- 
out efficiency correction for Ki obtained from signal 
MC. 

The uncertainty due to cross-feed of events between 
bins is conservatively estimated by taking the bias be- 
tween the fits with and without cross-feed. 

The uncertainty arising from the finite sample of flavor- 
tagged D — > Kgir + ir~ decays is evaluated by varying the 
numbers of flavor-tagged events in bins Ki within their 
statistical errors. 

The final results for x,y are corrected for the fit bias 
obtained from the fits of MC pseudo-experiments. The 
uncertainty due to the fit bias is taken from the difference 
of biases for various input values of x and y. 

The uncertainty due to errors of Ci and s» parameters is 
obtained by smearing the Ci and Sj values within their to- 
tal errors and repeating the fits for the same experimental 
data. We have performed a study of this procedure using 
the MC pseudo-experiments and analytical calculations. 
We find that the uncertainty obtained this way is sample- 
dependent for small B data samples and its average scales 
inversely proportional to the square root of sample size. 
It reaches a constant value for large B data samples (in 
the systematics-dominated case). This explains a some- 
what higher uncertainty compared to the CLEO estimate 
given in fl2| obtained in the limit of very large B sam- 
ple. In addition, the uncertainty in x, y is proportional 
to tb, and thus the uncertainty of the phases and 8b 
is independent of tb- As a result, the uncertainty of x, y 
in the S ± — > DK^ sample fit is 3-4 times larger than 
for B ± -> Dtt*. 



X. RESULTS FOR fo, r B AND 8 B 

We use frequentist treatment with the Fcldman- 
Cousins ordering to obtain the physical parameters 
H = (</>3,rB,^s) from the measured parameters z = 
[x- , y_, sg+, y+) as was done in the previous Belle analy- 
ses. In essence, the confidence level a for a set of physical 
parameters fi is calculated as 

ct(fi) = J p(z\/i)dz I J p(z\n)dz , (16) 

X>0) ' oo 

where p(z\n) is the probability density to obtain the mea- 
surement result z given the set of physics parameters \x. 
The integration domain T>(p) is given by the likelihood 
ratio (Fcldman-Cousins) ordering: 

p(z\nb cs t{z)) p(zo|/"bcst(2o))' 

where /J-best(z) is /i that maximizes p(z\fi) for the given 
z, and zq is the result of the data fit. 
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TABLE V. Systematic errors of x, y measurement for B — > D-k^ and B — > DK ± samples in units of 10 



Source of uncertainty 


B ± -i 






B ± ^ 


DK± 




Ax- 


Ay- 


Ax+ 


Ay+ 


Ax- 


Ay. 


Ax+ 


Ay+ 


Signal shape 


0.9 


1.9 


1.1 


5.0 


7.3 


7.4 


7.3 


5.1 


u, d, s, c continuum background 


0.9 


1.4 


0.8 


1.3 


6.7 


5.6 


6.6 


3.2 


BB background 


3.3 


1.6 


4.5 


1.1 


7.8 


12.2 


7.2 


6.1 


B — > Dir^ background 


— 


— 


— 


— 


1.2 


4.2 


1.9 


1.9 


Dalitz plot efficiency 


3.0 


1.9 


3.2 


1.6 


4.8 


2.0 


5.6 


2.1 


Cross-feed between bins 


0.4 


3.0 


0.7 


0.9 


0.4 


9.0 


0.0 


3.0 


Flavor-tagged statistics 


1.7 


2.0 


1.6 


2.0 


1.5 


2.7 


1.7 


1.9 


Fit bias 


0.4 


0.5 


0.4 


0.5 


3.2 


5.8 


3.2 


5.8 


Cj, Si precision 


2.6 


6.5 


2.6 


6.5 


10.1 


22.5 


7.2 


17.4 


Total without Ci,Si precision 


5.0 


5.0 


6.0 


6.0 


14.0 


19.4 


14.0 


11.3 


Total 


5.6 


8.2 


6.5 


8.8 


17.3 


29.7 


15.7 


20.7 
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<J> 3 (degrees) fa (degrees) 

FIG. 9. Two-dimensional projections of confidence region 
onto (03,<5b) and (03, rs) planes (one-, two-, and three stan- 
dard deviations). 



The difference with the previous Belle analyses is that 
the probability density p(z\fi) is a multivariate Gaussian 
PDF with the errors and correlations between x± and y± 
taken from the data fit result. In the previous analyses, 
this PDF was taken from MC pseudo-experiments. 

As a result of this procedure, we obtain the confidence 
levels (CL) for the set of physical parameters 03, r_s, 63- 
The confidence levels for one and two standard deviations 
are taken at 20% and 74% (the case of three-dimensional 
Gaussian distribution). The projections of the 3D sur- 
faces bounding one and two standard deviations volumes 
onto 03 variable, and (03, re) and (03, #g) planes are 
shown in Fig. |H1 

Systematic errors in \i are obtained by varying the mea- 
sured parameters z within their systematic errors (Gaus- 
sian distribution is taken) and calculating the RMS of 
Mbestt^)- In this calculation we assume that the system- 
atic errors are uncorrelated. In the case of Cj, Sj system- 
atica, we test that assumption: when the fluctuation in a 
and Si is generated, we perform the fits to both B + and 
B~ data with the same fluctuated Cj,Sj. We observe no 



significant correlation between resulting X- and x+ (y_ 
and y+). 

The final results are: 

03 = (77.31^ ±4.2 ±4.3)° 

r B =0.145 ±0.030 ±0.011 ±0.011 (18) 

8 B = (129.9 ± 15.0 ± 3.9 ± 4.7)°, 

where the first error is statistical, the second is systematic 
error without Cj , Sj uncertainty, and the third error is due 
to Ci , Si uncertainty. 

We do not calculate the statistical significance of CP 
violation as it is done in the previous analyses by tak- 
ing the CL for 03 = 0: this number is purely based on 
the behavior of the tails of p(z\fi) distribution far from 
the central value, and Gaussian assumption can lead to 
overestimation of CP violation significance. As a prelim- 
inary number we use the estimate of probability of the 
fluctuation in the difference of number of events in bins 
for B + and B~ data: the probability of such fluctuation 
in the case of CP conservation is p = 0.42%. 

XI. CONCLUSION 

We report the results of a measurement of the unitarity 
triangle angle 03 using a model-independent Dalitz plot 
analysis of D — >• 7\^7r + 7r _ decay in the process B ± —> 
DK . The measurement was performed with a full data 
sample of 711 ftr 1 (772 x 10 6 BB pairs) collected by 
the Belle detector at T(45). The model independence 
is reached by binning the Dalitz plot of D — > Kgir + ir~ 
decay and using the strong phase coefficients for bins 
measured by CLEO experiment [l2[ . We obtain the value 
</> 3 = (77.3± 1 |;g±4.2±4.3)°; of the two possible solutions 
we choose the one with < 03 < 180°. We also obtain 
the value of the amplitude ratio r B = 0.145 ± 0.030 ± 
0.011 ±0.011. These results are preliminary. 
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This analysis is a first application of the novel method 
of 03 measurement. Although currently it does not offer 
significant advantages over the model-dependent Dalitz 
plot analyses of the same decay chain, it is promising for 
the measurement at super-B factories [l8|,[l!|. We expect 
that the statistical error of the 03 measurement using the 
statistics of 50 ab _1 to be available at the super-B fac- 
tory will reach 1-2°. With the use of BES-III data 
the error due to the phase terms of D — > KgTr + ir~ decay 
will decrease to 1° or less. We also expect that the exper- 
imental systematic error can be kept at the level below 
1° since most of its sources are limited by the statistics 
of the control channels. 
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